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In view of the significantly different proton charge radius extracted from muonic hydrogen Lamb shift mea¬ 
surements as compared to electronic hydrogen spectroscopy or electron scattering experiments, we study in this 
work the photoproduction of a lepton pair on a proton target in the limit of very small momentum transfer as a 
way to provide a test of the lepton universality when extracting the proton charge form factor. By detecting the 
recoiling proton in the -yp —> l~l + p reaction, we show that a measurement of a ratio of e _ e + + pi~ pi + over 
e“e + cross sections with an absolute precision of 7 x 1CL 4 , would allow for a test to distinguish, at the 30- 
level, between the two different proton charge radii currently extracted from muonic and electronic observables. 


Recent extractions of the proton charge radius Re from 
muonic hydrogen Lamb shift measurements m m are in 
strong contradiction, by around 7 standard deviations, with 
the values obtained from energy level shifts in electronic hy¬ 
drogen |0 or from electron-proton elastic scattering |j4)|5l : 

HP O : R e = 0.8409 (4) fm, (1) 

ep 13 : R e = 0.8775 (51) fm. (2) 

This so-called ’’proton radius puzzle” has triggered a large ac¬ 
tivity and is the subject of intense debate, see e.g. BH for 
recent reviews, and references therein. Lepton universality re¬ 
quires the same radius to enter the electronic and muonic ob¬ 
servables. If the different Re extractions cannot be explained 
by overlooked corrections, it would point to a violation of 
electron-muon universality. Several scenarios of new, beyond 
the Standard Model, physics have been proposed by invoking 
new particles which couple to muons and protons, but much 
weaker to electrons, see e.g. mm. Such models would also 
lead to large loop corrections to the muon’s anomalous mag¬ 
netic moment, (g — 2) /( , which presently displays a 3<r devi¬ 
ation between experiment and its Standard Model prediction, 
see e.g. m. Explaining both the (g — 2) ; , discrepancy and 
the proton radius puzzle by new particles coupling mainly to 
muons seems an attractive perspective. It does however re¬ 
quire a significant fine-tuning, especially for larger values of 
the conjectured new particle masses J3- To test the electron- 
muon universality, it has been proposed by the MUon proton 
Scattering Experiment (MUSE) 03 to make a simultaneous 
measurement of both pip and ep elastic scattering, extracting 
the proton charge form factor from both measurements. Be¬ 
sides the plans to measure pip elastic scattering, several new 
experiments are underway to extend ep scattering to lower 
momentum transfer values, down to 10 1 GeV 2 , and to cross¬ 
check its systematics naira. All of these tests require abso¬ 
lute cross section measurements, with a required precision on 
each absolute cross section at the level of 1 % or better. 

In order to reach in pi scattering experiments the precision 
on Re comparable with e" scattering experiments, indicated 
in Eq. Q. we propose in this work a new experimental avenue 
through a relative cross section measurement of the photo¬ 
production of e _ e + versus pi~ /x + pairs on a proton target. Be¬ 
sides being a well studied process El, the photo-production 
of a lepton pair has the advantage that one produces e _ e + and 
pi~ pi + final states with the same beam, and thus the overall 


normalization uncertainty drops out of their ratio. The analy¬ 
sis presented in this work shows that through a detection of the 
recoiling proton in the 7 p -A l~l + p reaction, a measurement 
of the ratio of the e _ e + cross section below pi~ pi + thresh¬ 
old versus the e _ e + + pi~ pi + cross section sum above pi~ pi + 
threshold with an absolute precision of around 7 x 1CT 4 will 
allow to distinguish, at the 3er level, between the proton Re 
extractions from muonic and electronic observables. Further¬ 
more, we show that the linear photon polarization asymme¬ 
try has a discriminatory power between the e _ e + and pi~ pi + 
channels in the pi~ pi + threshold region. 




FIG. 1: Bethe-Heitler mechanism to the 7 p —> l~l + p process, where 
the four-momenta of the external particles are: k for the photon, 
p(p') for initial (final) protons, and Z_, l + for the lepton pair. 

We will consider the lepton pair production on a proton 
target, 7 p -A l~l + p, in the limit of very small momentum 
transfer, defined as t = (p — p') 2 , with four-momenta as in¬ 
dicated on Fig. |T] Furthermore, we will use in the following 
the Mandelstam invariant s = ( k + p ) 2 = M 2 + 2 ME 1 , with 
M the proton mass and E 1 the photon lab energy, as well 
as the squared invariant mass of the lepton pair, defined as 
M ; 2 = (Z_ + l+) 2 . In the limit of small — t, the Bethe-Heitler 
(BH) mechanism, shown in Fig. [T| to tally dominates the cross 
section of the 7 p -A l~l + p reaction. As the momentum trans¬ 
fer t is the argument appearing in the form factor (FF) in the 
BH process, a measurement of the cross section in this kine¬ 
matic regime will allow to access the proton electric FF Ge p 
at small spacelike momentum transfer, with a very small con¬ 
tribution of the proton magnetic FF G m p . As the differential 
cross section for the BH process is strongly peaked for lep¬ 
tons emitted in the incoming photon direction, and as we aim 
to maximize the BH contribution in this work in order to ac¬ 
cess Ge p , we will study the 7 p -A l~l + p process when (only) 
detecting the recoiling proton’s momentum and angle, thus ef¬ 
fectively integrating over the large lepton peak regions. The 
lab momentum of the proton is in one-to-one relation with 
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the momentum transfer t: \p'\ lab = 2My/r(l + r), with 
r = — t/{AM 2 ). Furthermore, for a fixed value of t, the 
recoiling proton lab angle 0 ^ is expressed in terms of in¬ 
variants as : 


cos 0 


lab 

P 


M 2 + 2(s + M 2 )t 
2 (s - M 2 )i/r(l 4 - r) 


(3) 


with a = e 2 /At: ps 1/137, where /3 = ^/l — is the 

lepton velocity in the Z“Z + c.m. frame, with m the lepton 
mass, and where the proton FFs Ge p and Gm p are functions 
of t. The weighting coefficients multiplying the FFs in Eq. ([4]) 
have the following general structure : 


The differential cross section for the dominating BH pro¬ 
cess to the 7 p —> l~l + p reaction has been studied in different 
contexts in the literature Ifl8l - l20l . In this work, we will con¬ 
sider the cross section differential in the momentum transfer 
t and invariant mass of the lepton pair M'/j , and integrated 
over the lepton angles, which corresponds with detecting the 
recoiling proton only. This cross section can be written as : 

da BH a 3 4/3 1 

dtdM 2 = {s-M 2 ) 2 ' t 2 (M 2 - t ) 4 ' 1 T t 

x | Ce G 2 Ep + Cm t G 2 Mp } , (4) 


Ce,m = C$ M + In , (5) 


where the second term expresses the large logarithmic en¬ 
hancement in the limit of small lepton mass in the BH process. 
The coefficients C^ M , and C E \ M are found to be expressed 
through invariants as : 


4 1} 

= t (s - 

- M 2 ) (s-M 2 - Ml 

+ t) [Ml + 6Mlt + t 2 + Am 

2 mI\ 



+ (M 2 

- t ) 2 [t 2 Ml + M 2 (Ml + t ) 2 + 4to 

2 m 2 m 2 \ , 


(6) 


= 

s — M 2 ) (s-M 2 -M 2 +t) [M& + 

t 2 + Am 2 (Ml 

+ 2 1 — 2m 2 )] 



+ (M 2 

— t) 2 [—M 2 (Ml +1 2 

) + 2m 2 (— t 2 

- 2 M 2 Ml + 4 m 2 M 2 )] , 

(7) 

/'f(l) 

C 'M 

= 4 1} 

-2M 2 (1+t) (M 2 - 

t) 2 [M/ l + t 2 

+ 4 m 2 Ml\ , 


(8) 
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+ 2M 2 (1 +t) (Ml - 

t) 2 [Mp + t 2 

+ 4m 2 (M; 2 — t 

— 2?77, 2 )] . 

(9) 


In Fig. El we show the differential cross section da/dt dM'/j 
for 7 p ~^{l~l + )p which is accessed by measuring the recoil¬ 
ing proton’s momentum and angle, for E 1 = 0.5 GeV and 
for three values of — t : —t = 0.01 GeV 2 (corresponding 
with recoil proton momentum \p'\ lab = 100 MeV/c), —t = 
0.02 GeV 2 (\p'\ lab = 142 MeV/c), and -f = 0.03 GeV 2 
(| p'\ lab = 174 MeV/c). As t is the argument entering the pro¬ 
ton FFs, the values shown are chosen to cover the lower range 
of the high-precision elastic ep scattering experiments ma, 
as well as the values for which the future MUSE elastic pp 
scattering experiment |fj~5l plans to take data. Furthermore, 
the differential cross section da/dtdM'/ is shown as func¬ 
tion of the squared lepton invariant mass M/, which is dialed 
through the recoiling proton lab angle, as shown on the lower 
panel of Fig. [ 2 ] We show the cross section in a range of M/, 
which is kinematically separated from background channels, 
well above the Compton and 7r° production processes on a 
proton, corresponding with sharp peaks at M/ = 0 and at 
M ; 2 = 0.018 GeV 2 respectively, and below the threshold for 
7 T 7 r production, which starts at M ; 2 = 0.078 GeV 2 . One no¬ 
tices from Fig. [ 2 ] that around M z 2 = 0.06 GeV 2 , the p~p + 
cross section is around a factor of 10 smaller than the e“e + 
cross section, and increases with increasing M/[. By measur¬ 
ing the cross section through detecting the recoiling proton 
momentum and angle in the M/ window above the 7r° peak 


and below irn threshold, and comparing cross sections at a 
fixed value of t above and below p~ p + thresholds, it opens 
the possibility for a high-precision extraction of the cross sec¬ 
tion ratio : 


R p /e — 


da(p p + + e e + ) 
da(e~e + ) 


GO) 


where da stands for da/dt dM 2 v The potential advantage of 
such a ratio measurement is that absolute normalization un¬ 
certainties to first approximation drop out. Indeed, at a fixed 
value of t, the e _ e + cross section can be fixed by measur¬ 
ing the cross section below p~ p + threshold, and the corre¬ 
sponding normalization, due to Ge p , can be used to deter¬ 
mine the e _ e + cross section above p~ p + threshold. A subse¬ 
quent measurement of the sum of e _ e + + p~ p + cross sec¬ 
tions above p~ p + threshold, then allows to extract the ra¬ 
tio R p / e , which is displayed in Fig. [3] One sees that in the 
kinematic range where only the e~er and p~p + channels 
are contributing, this ratio varies between 1.10 to 1.14. We 
like to notice that corrections, notably radiative corrections, 
to first order also drop out of this ratio, measured at the same 
value of the recoiling proton momentum and angle. An ac¬ 
curate measurement of this ratio can therefore be envisaged, 
opening a new perspective to perform a test of lepton univer¬ 
sality. We have demonstrated this sensitivity in Fig. [3] by 
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FIG. 2: Upper panel: comparison of the (lepton pair) invariant mass 
dependence of the 7 p —> e~e + p process (upper three curves) vs the 
7 p —> p~ p + p process (lower three curves) at E 1 = 0.5 GeV, and 
for three values of the momentum transfer t as indicated. The lower 
panel shows the corresponding kinematic relation between the lepton 
pair invariant mass and the recoiling proton lab angle. 


varying the electric FF value entering the e e + production 
process, denoted by G e Ep , which we take from 0151, versus 
the electric FF value entering the p + production process, 
denoted by G Ep . We compare the case of lepton universal¬ 
ity G Ep = G Ep , with a lepton universality violation scenario 
in which G Ep /G e Ep = 1.01. The latter value is motivated 
by the around 1% larger value of G Ep at —t = 0.03 GeV 2 , 
resulting from the proton radius as extracted from the muonic 
hydrogen Lamb shift of Eq. 0 - in comparison with the proton 
radius as extracted from ep scattering or electronic hydrogen 
Lamb shift, given by Eq. 0 - We like to note that more recent 
studies, using a more conservative treatment of systematic er¬ 
rors in elastic ep scattering, extract a radius in agreement with 
Eq. Q, but with twice larger error GTj . 

We see from Fig.[3]that these two scenarios for the proton 
radius lead to a difference in the cross section ratio R p / e of 
Eq. (10 1 of around 2 x 10 -3 . A measurement of this ratio 
R p / e with a precision of around 7 x 10 ” 4 would therefore 
provide a test of such difference at the 3cr level. In Fig. [3] 
we show the 1,3, and 5 a error bands corresponding with a 
a of 7 x 10 -4 . Although the Compton and 7r° photoproduc¬ 
tion backgrounds have been eliminated by our choice of the 
kinematical range in , we also estimated the background 
contamination due to the 7 p -7 ( 7 r° 7 )p process l 22 l in the 
relevant range M/) ss 0.07 GeV 2 . We found that in order for 


the contamination on R p / e to be bounded by AR p / e < 10 -3 , 
the required precision on the 7 p —> 7 r° 7 p cross section should 
be in the 10 - 30 % range, which has already been achieved 
by present experiments. Furthermore, the remaining concern 
is the physical background due to the indistinguishable time¬ 
like Compton scattering process, which results in the same fi¬ 
nal state. We estimate this timelike Compton process at the 
relatively low energy and momenta considered here by its 
Born contribution, corresponding to a nucleon intermediate 
state l23l . The inclusion of the timelike Compton Born con¬ 
tribution is also shown in Fig. [3] We found the effect due to 
the Born contribution in most kinematics around a factor of 5 
smaller than the effect due to the 1 % variation in the value 
of G Ep , shown in Fig. |3j Although the timelike Compton 
process requires further study, we like to notice that it is not 
totally unknown, and a moderate accuracy in the 30 % range 
would be sufficient for its contamination to be bounded by 
AR p / e < 10 -3 . We are therefore confident that the ratio of 
Eq. ( [Tol l i s a promising observable for a lepton universality 
test at the 1 % level between G Ep and G Ep . Besides provid¬ 
ing such a test, a measurement of the absolute 7 p -7 (e~e + )p 
cross section below p~p + threshold will provide a useful 
cross-check on the extraction of Ge p from elastic ep scatter¬ 
ing. As the systematics in both reactions are quite different, 
an independent measurement of Ge p may yield a valuable fur¬ 
ther clue to shed light on the “proton radius puzzle”. 



/W,| (GeV 2 ) 

FIG. 3: Ratio R p / e of the 7 p — > (e“e + + p~p + )p vs 7 p — >■ 
(e“e + )p differential cross sections, according to Eq. ( | 10| ). The 
lower (blue) curve corresponds with the lepton universality result, 
for which G Ep = G Ep . The upper (red) curve assumes a violation 
of lepton universality, for which G Ep /G e Ep = 1.01. The blue bands 
denote 1, 3, 5 a bands corresponding with a a of 7 x 10~ 4 around 
the blue curve. The difference between the blue curve and the dotted 
curve is an estimate of the physical background due to the contribu¬ 
tion of the timelike Compton process. 

Besides the unpolarized cross section , we may also con¬ 
sider the sensitivity of polarization observables to distinguish 
between the e _ e + and p~p + production processes. We will 
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study here the case of the linear photon asymmetry defined as: 


dcr|| — dcr± 
dan + da± ’ 


( 11 ) 


where da\\ ( da± ) stands for the differential cross section for a 
photon with linear polarization parallel (perpendicular) to the 
plane spanned by the photon and recoiling proton momenta. 
When measuring the recoiling proton only, the asymmetry 
above //T //+ threshold is given by the following weighted 
sum of the asymmetries of the e _ e + and /i + channels: 


Aiu(e e + + ii n + ) = 


R 


'M/e 


{Aiu(e e + ) 


+ ( i? M/e - 1 ) A iu (ii-il+)} . ( 12 ) 

We show in Fig. [4] the linear photon asymmetry in the kine¬ 
matic range around n~ n + threshold. It is seen that the linear 
photon asymmetry is very small for the e _ e + channel. How¬ 
ever, for the h~il + channel the asymmetry reaches a value 
approaching — 1 at threshold and decreases in absolute 

value by going away from the threshold. Such behavior arises 
due to an exact cancellation, at threshold, between the 

analytical and lepton mass logarithmic terms in the analogous 
expression as Eq. (|5j» for the contribution proportional to G 2 Ep 
to a | . In the whole Mfi range of interest, the asymmetry for 
fi~ IL + production takes on large values as can be seen from 
Fig. 0 A direct measurement of the // /i + asymmetry may 
therefore give a clear tool to separate the two channels. If 
the lepton pair is undetected, and only the recoiling proton is 
measured, the measured asymmetry above threshold is 

diluted by the fi~ fi + /e~e + ratio, as given by Eq. 
reaches values around -5 % as can be seen from Fig. 


12 | 

|4T 


and 



FIG. 4: Linear photon asymmetry Aiu of the 7 p —» ( l~l + )p pro¬ 
cess. The dashed (blue) curve corresponds with e“e + production; 
the dashed-dotted (red) curve corresponds with it~ ii + production. 
The solid (black) curve is the asymmetry corresponding with the sum 
of the e“e + + /r _ M + channels according to Eq. 1 12 1 . 


In conclusion, in view of the sizably different proton charge 
radius presently extracted from electronic and muonic observ¬ 
ables, we proposed a new lepton universality test which meets 


the required precision goal to distinguish between them. We 
demonstrated that such a test is possible by comparing the 
photoproduction of a lepton pair on a proton, through detec¬ 
tion of the recoiling proton. We showed that the measurement 
of the ratio of photoproduction cross sections of e _ e + +/r _ /r + 
vs e _ e + pairs with an absolute precision of around 7 x 10 1 
will allow to distinguish, at the 3<r level between the differ¬ 
ent proton Re extractions from muonic and electronic ob¬ 
servables. Such an experiment can be performed at existing 
electron facilities such as the Mainz Mikrotron (MAMI) and 
Jefferson Lab, thus adding a further piece of evidence towards 
an understanding of the “proton radius puzzle”. 
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